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(57) Abstract: Signal recovery in functional magnetic resonance imaging (fMRI) is provided by generating a single excitation pulse 
£j and exciting a target region of a subject with the generated excitation pulse. A first image is obtained using a first partial k-space 

frame of the target region. A compensation pulse is generated and the target region excited with the compensation pulse. A second, 

compensated, image is obtained subsequent to the excitation by the compensation pulse using a second partial k-space frame of 
£^ the target region. The first and second images are combined to form a combined image of the target region. The first and second 

obtaining steps are carried out sequentially during a single quadratic excitation pulse. 
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SINGLE-SHOT EPI WITH SIGNAL RECOVERY FROM SUSCEPTIBILITY-INDUCED LOSSES 

Claim of Priority 

5 The present application claims priority from United States Provisional Patent 

Application Serial No. 60/241,999 filed October 20, 2000, the disclosure of which is 
incorporated by reference herein as if set forth fully herein 

Field of the Invention 

The present invention is related to functional magnetic resonance imaging 
10 (fMRI) methods. 

Background of the Invention 

Functional magnetic resonance imaging (fMRI) has been used widely in brain 
imaging studies for the past several years. See e.g. , J. W. Belliveau, D. N. Kennedy, 
R. C. McKinstry, B. R. Buchbinder, R. M. Weisskoff, M. S. Cohen, J. M. Vevea, T. J. 

15 Brady and B. R. Rosen, "Functional mapping of the human visual cortex by magnetic 
resonance imaging," Science 254, 716-719, 1991; K. K. Kwong, J. W. Belliveau, D. 
A. Chesler, L E. Goldberg, R. M. Weisskoff, B. P. Poncelet, D. N. Kennedy, B. E. 
Hoppel, M. S. Cohen, R. Turner, H.-M. Cheng, T. J. Brady and B. R. Rosen, 
"Dynamic magnetic resonance imaging of human brain activity during primary 

20 sensory stimulation," Proc. Natl Acad. Sci. USA 89, 5675-5679, 1992; P. A. 

Bandettini, E. C. Wong, R. S. Hinks, R. S. Tikofsky, and J. S. Hyde, "Time course 
EPI of human brain function during task activation," Magn. Resort. Med. 25, 390-397, 
1992; S. Ogawa, D. W. Tank, R. Menon, J. M. Ellerman, S.-G. Kim, H. Merkle, and 
K. Ugurbil, "Intrinsic signal changes accompanying sensory stimulation: functional 

25 brain mapping with magnetic resonance imaging," Proc. Natl Acad. Set USA 89, 
5951-5955, 1992; and R. Menon, S. Ogawa, D. W. Tank, and K. Ugurbil, "4 Tesla 
gradient recalled echo charateristics of photic stimulation-induced signal changes in 
the human primary visual cortex," Magn. Resort. Med. 30, 380-386, 1993. One of the 
most often used methods is the gradient-recalled echo-planar imaging (EPI) technique 

30 because of its good sensitivity to the blood oxygenation level dependent signal and 

high speed. See S. Ogawa, R. S. Menon, D. W. Tank, S.-G. Kim, H. Merkle, J. M. Ellerman 
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and K. Ugurbil, "Functional brain mapping by blood oxygenation level-dependent contrast 
magnetic resonance imaging," Biophys. J. 64, 803-812, 1993. However, its usage is 
limited in areas with severe static inhomogeneity induced by susceptibility effect near 
air/tissue interfaces. 

5 One potential problem in gradient-recalled EPI vising a long echo time is the 

severe signal losses at areas with large static inhomogeneities. These areas include 
the ventral frontal, medial temporal and inferior temporal regions that experience 
inhomogeneities induced by the susceptibility effects near air/tissue interfaces. For 
fMRI studies that use both the gradient-recalled EPI and high field scanners, these 
10 signal losses may prevent investigation of the human cognitive processes such as the 
memory and attention studies. Methods have been developed to recover the signal 
losses, however, these methods typically involved multiple excitations, thus, 
compromising the temporal resolution. 

Susceptibility artifacts can be manifested primarily in two ways: signal losses 

1 5 and geometric distortions. In general, a long echo time makes an MRI system more 
prone to signal losses in the presence of an inhomogeneous field because of the intra- 
voxel dephasing, and a long readout time typically leads to geometric distortions due 
to the reduced sampling frequency and reduced readout gradient strength. Typically, 
pronounced susceptibility-related field variation along the slice-selective direction in 

20 combination with a long echo time results in severe signal losses, while the 

inhomogeneity in-plane combined with a long readout time leads to geometric 
distortions. Thus, signal losses can be caused by the susceptibility-induced gradient 
along the slice-selective direction. Because of the long echo time typically used in 
fMRI experiments to maximize the sensitivity toward the signal changes, the signal 

25 losses at areas near air/tissue interfaces may be severe. Refined methods to recover 
these signals may be needed in order to study brain function at these areas. 

Several research groups have addressed these sorts of signal losses using 
various techniques. One such technique is to use a thinner slice thickness to reduce 
the field change across the slice. See L R. Young, I. J. Cox, D. J. Bryant, and G. M. 

30 Bydder, "The benefits of increasing spatial resolution as a means of reducing artifacts 
due to field inhomogeneities," Magn Resoru Imag. 6, 585-590, 1988. This technique 
may be relatively easy to implement but it may reduce SNR as well as the spatial 
coverage per unit time. 
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Frahm et al. originally proposed to use multiple refocusing gradients to 
effectively compensate the field inhomogeneities. J. Frahm, K. D. Merboldt, W. 
Hanicke, "Direct FLASH MR imaging of magnetic field inhomogeneities by gradient 
compensation," Magn. Reson. Med. 6, 474-480, 1988. This method was later adopted 
5 by several other groups and applied more recently in functional MRI. See e.g., R. J. 
Ordidge, J. M. Gorell, J. C. Deniau, R. A. Knight, J. A. Helpern, "Assessment of 
relative brain iron concentrations using T2-weighted and T2*-weighted MRI at 3 
Tesla," Magn. Reson. Med. 32, 335-341, 1994; R. T. Constable, "Functional MR 
imaging using gradient-echo echo-planar imaging in the presence of large static field 

10 inhomogeneities," J. Magn. Reson. Im. 5, 746-752, 1995; Q. X. Yang, B. J. 

Dardzinski, S ? Li, P. J. Eslinger, M. B. Smith, "Multi-gradient echo with susceptibility 
compensation (MGESIC): demonstration of fMRI in the olfactory cortex at 3T," 
Magn. Reson. Med. 37, 331-335, 1997; R. T. Constable, D. D. Spencer, "Composite 
image formation in z-shimmed functional MR imaging," Magn. Reson. Med 42, 1 1 0^ 

15 117, 1999; and V. A. Stenger, F. E. Boada, and D. C. Noll, "Gradient compensation . 
method for the reduction of susceptibility artifacts for spiral fMRI data acquisition," 
Proc. ISMRM, p. 538, 1999. 

Because the superimposed gradient field across the slice is often not linear, 
one compensatory gradient is generally not sufficient to compensate the entire slice. 

20 Thus, a set of linear gradients is typically needed to compensate the nonlinear field 
segment-by-segment to achieve satisfactory results. When the number of the linear 
gradients increase, i.e., increments become finer, the nonlinear field can be better 
compensated. Despite the effectiveness in recovering signal, the time-consuming 
nature of such techniques may limit their practical value in routine fMRI experiments. 

25 In practice, as many as sixteen repetitions may be needed to sum up to a uniform 
image. Most of cognitive fMRI experiments cannot be performed this way. 

More sophisticated methods were also proposed that showed promise in 
reducing the number of compensating gradients to a much more tolerable level. The 
efficiency is much increased by using high-order field compensation. Cho et al. 

30 proposed tailored pulse with a quadratic profile that has shown improved tolerance 
toward field inhomogeneity. Z. H. Cho, and Y. M. Ro, "Reduction of susceptibility 
artifact in gradient-echo imaging," Magn. Reson. Med. 23, 193-196, 1992. Glover et 
al also presented a method using high order phase compensation by obtaining a field 
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profile for each subject and incorporating it into the phase profile of the excitation 
pulse. G. Glover, S. Lai, "Reduction of susceptibility effects in fMRI using tailored 
RF pulses," Proc. ISMRM, p.298, 1998. A similar concept was used in a recent report 
using a two-shot technique with explicitly matched RF excitation. N. K. Chen, A. M. 
5 Wyrwicz, "Removal of intravoxel dephasing artifact in gradient-echo images losing a 
field-map based RFrefocusing technique," Magn. Reson. Med. 42, 807-812, 1999. 
Another recent method used a two-shot technique combining a quadratic excitation 
pulse and the compensatory gradient. J. Mao, and A. W. Song, "Intravoxel rephasing 
of spins dephased by susceptibility effect for EPI sequences," Proc. ISMRM, p. 1982, 

10 1 999. The resultant phase profile can be used to better match the susceptibility- 
induced gradients when an appropriate compensatory gradient is used. The two 
excitations can be implemented back-to-back within one run to allow fMRI 
experiments to be carried out; however, the effective repetition time is still doubled. 
Images from the two excitations can then be combined to achieve uniform spatial 

15 coverage across the inhomogeneous areas. 

Summary of the Invention 

The present invention provides systems, methods and computer program 
20 products which provide fMRI signal recovery from a single excitation. Such 

methods, systems and computer program products may be particularly suitable and 
useful for fMRI studies and applications carried out in or about inhomogeneous areas 
with high temporal resolution. 

In particular embodiments of the present invention, signal recovery in 
25 functional magnetic resonance imaging (fMRI) is provided by generating a single 

excitation pulse and exciting a target region of a subject with the generated excitation 
pulse. A first image is obtained using a first partial k-space frame of the target region. 
A compensation pulse is generated and the target region subjected to the 
compensation pulse. A second, compensated, image is obtained subsequent to the 
30 target region being subjected to the compensation pulse using a second partial k-space 
frame of the target region. The first and second images are combined to form a 
combined image of the target region. The first and second obtaining steps are carried 
out sequentially dining a single quadratic excitation pulse. 
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In particular embodiments of the present invention, the single excitation pulse 
is a matched quadratic excitation pulse. Furthermore, the compensation pulse may be 
a z-shimming pulse. 

In further embodiments of the present invention, the first partial k-space frame 
5 has an associated first sampling direction and the second partial k-space frame has an 
associated second sampling direction. Preferably, the first sampling direction and the 
second sampling direction are substantially the same direction. 

In still additional embodiments of the present invention, the first partial k- 
space frame and the second partial k-space frame are sampled so that a center of each 
10 partial k-space frame is proximate. In such embodiment, the first image may be 
obtained by completely sampling a first half of the first partial k-space frame and 
partially sampling a second half of the first partial k-space frame. The second image 
is then obtained by partially sampling a first half of the second k-space frame and 
completely sampling a second half of the second k-space frame. Preferably, the 
15 sampling of the first half and the second half of the first k-space frame and the 

sampling of the first half and the second half of the second k-space frame are carried 
out in substantially the same direction. 

In additional embodiments of the present invention, the Echo Time (i.e., the 
time difference between the excitation pulse and the center of k-space acquisition) 
20 (TE) associated with the first partial k-space frame and the TE associated with the ..• 
second partial k-space frame are asymmetric with reference to the beginning of the 
respective partial k-space frames and are mirrored about a division between the first 
and second partial k-space frames. 

In still further embodiments of the present invention, alignment indicia are 
25 inserted in the first and second images and the alignment indicia are aligned in each of 
the first and second images when the images are combined. In particular, the 
alignment indicia may be a centerline of the first k-space frame and a centerline of the 
second k-space frame. 

In yet other embodiments of the present invention, a centerline of the first 
30 partial k-space frame is sampled in two sampling directions and a centerline of the 
second partial-k-space frame is also sampled in two sampling directions. 

In additional embodiments of the present invention, the target region of a 
subject is an in vivo ventral frontal or inferior temporal area of the human brain. 
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Furthermore, the combined image may provide information about the function of the 
human brain including human memory and attention processes in areas near air/tissue 
interfaces. 

While embodiments of the present invention are described above with 
5 reference to methods, as will be appreciated by those of skill in the art, embodiments 
of the present invention may also be provided as systems and/or computer program 
products. 

Brief Description of the Drawings 
10 Figure 1 is a block diagram of an fMRI system according to embodiments of 

the present invention. 

Figure 2A is a flowchart illustrating operations according to embodiments of 
the present invention. 

Figure 2B is a flowchart illustrating operations according to further 
1 5 embodiments of the present invention. 

Figure 3A is a schematic illustration of a single-shot EPI sequence according 
to embodiments of the present invention. 

Figure 3B is a schematic illustration of a resultant k-space trajectory using the 
sequence shown in Figure 3 A according to embodiments of the present invention. 
20 Figures 4A and 4B are screen images of single-shot axial echo-planar images 

of a water phantom obtained according to methods of the present invention. Figure 
4A corresponds to images reconstructed from a first partial k-space and Figure 4B 
corresponds to images reconstructed from the second partial k-space. 

Figures 5A-5C are screen images of axial images of a human brain using EPI 
25 acquisitions. Figure 5A illustrates images obtained using conventional single shot 
gradient recalled echo planar images. Figure 5B illustrates images obtained using a 
two-shot acquisition with matched quadratic excitation and z-shim gradient. Figure 
5C illustrates images obtained using embodiments of the present invention utilizing a 
the single-shot acquisition with the same quadratic excitation and z-shim gradient. 

30 

Description of Embodiments of the Invention 

The present invention now will be described more fully hereinafter with 
reference to the accompanying drawings, in which illustrative embodiments of the 

6 
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invention are shown. This invention may, however, be embodied in many different 
forms and should not be construed as limited to the embodiments set forth herein; 
rather, these embodiments are provided so that this disclosure will be thorough and 
complete, and will fully convey the scope of the invention to those skilled in the art. 
5 As will be appreciated by one of skill in the art, the present invention may be 

embodied as a method, system, or computer program product. Accordingly, the 
present invention may take the form of an entirely hardware embodiment, an entirely 
software embodiment or an embodiment combining software and hardware aspects all 
generally referred to herein as a "circuit" or "module." Furthermore, the present 

1 0 invention may take the form of a computer program product on a computer-usable 
storage medium having computer-usable program code embodied in the medium. 
Any suitable computer readable medium may be utilized including hard disks, CD- 
ROMs, optical storage devices, a transmission media such as those supporting the 
Internet or an intranet, or magnetic storage devices. 

15 Computer program code for carrying out operations of the present invention 

may be written in an object oriented programming language such as Java®, Smalltalk 
or C++. However, the computer program code for carrying out operations of the 
present invention may also be written in conventional procedural programming 
languages, such as the "C" programming language. The program code may execute 

20 entirely on the user's computer, partly on the user's computer, as a stand-alone 

software package, partly on the user's computer and partly on a remote computer or 
entirely on the remote computer. In the latter scenario, the remote computer may be 
connected to the user's computer through a local area network (LAN) or a wide area 
network (WAN), or the connection may be made to an external computer (for 

25 example, through the Internet using an Internet Service Provider). Furthermore, the 
user's computer, the remote computer, or both, may be integrated into other systems, 
such as an MRI system. 

The present invention is described below with reference to flowchart 
illustrations and/or block diagrams of methods, apparatus (systems) and computer 

30 program products according to embodiments of the invention. It will be understood 
that each block of the flowchart illustrations and/or block diagrams, and combinations 
of blocks in the flowchart illustrations and/or block diagrams, can be implemented by 
computer program instructions. These computer program instructions may be 
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provided to a processor of a general purpose computer, special purpose computer, or 
other programmable data processing apparatus to produce a machine, such that the 
instructions, which execute via the processor of the computer or other programmable 
data processing apparatus, create means for implementing the functions/acts specified 
5 in the flowchart and/or block diagram block or blocks. 

These computer program instructions may also be stored in a computer- 
readable memory that can direct a computer or other programmable data processing 
apparatus to function in a particular manner, such that the instructions stored in the 
computer-readable memory produce an article of manufacture including instruction 

10 means which implement the function/act specified in the flowchart and/or block 
diagram block or blocks. 

The computer program instructions may also be loaded onto a computer or 
other programmable data processing apparatus to cause a series of operational steps to 
be performed on the computer or other programmable apparatus to produce a 

1 5 computer implemented process such that the instructions which execute on the 
computer or other programmable apparatus provide steps for implementing the 
functions/acts specified in the flowchart and/or block diagram block or blocks. 

Embodiments of the present invention provide for obtaining multiple 
acquisitions of MRI information from a single excitation pulse where a subsequent 

20 acquisition is compensated and the multiple acquisitions combined to provide an 
image. Although the conventional two-shot technique using matched quadratic 
excitation can be sufficient in achieving uniform images at inhomogeneous areas, the 
time penalty of such an approach may limit applications where both high temporal 
resolution (such as in event-related studies) and large spatial coverage are desired. 

25 Thus, embodiments of the present invention provide combining two acquisitions 
within one quadratic excitation which may improve the temporal resolution. 

One approach to combining two acquisitions from a single excitation is to 
sequentially acquire fc-space data twice within one excitation, with the second pass in 
£-space compensated with appropriate z-shim gradient. Two separate images can be 

30 reconstructed and combined to obtain a uniform spatial coverage. A drawback of this 
approach may be that a prolonged data acquisition window which may exceed the 
T2* value of the gray matter. Two unwanted features can result from this drawback: 
1) the echo times of the two acquisition are much different causing severe mixing of 
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the T2* contrast; and 2) the second acquisition may have very small signal intensity 
from the T2* decay, thus reducing the effectiveness of signal compensation. 

As described in more detail below, embodiments of the present invention may 
overcome such potential limitations through the novel use of partial &-space 
5 acquisition. Partial £-space acquisition is described in United States Patent No. 
4,767,991 entitled "METHOD OF HIGH SPEED IMAGING WITH IMPROVED 
SPATIAL RESOLUTION USING PARTIAL K-SPACE ACQUISITIONS* which issued 
August 30, 1988, the disclosure of which is incorporated herein by reference as if set 
forth fully herein. Partial £-space acquisition is also described in J. R. MacFall, N. J. . 

10 Pelc, and R. M. Vavrek, "Correction of spatially dependent phase shifts for partial 
Fourier imaging," Magn. Reson. Imaging. 6, 143-155, 1988 and A. Jesmanowicz, P. 
A. Bandettini, J. S. Hyde, "Single-shot half k-space high-resolution gradient-recalled 
EPI for fMRI at 3 Tesla," Magn. Reson. Med 40, 754-762, 1998. 

A system incorporating embodiments of the present invention is illustrated in 

15 Figure 1. As seen in Figure 1, an fMRI system 10 includes an MRI control system- 
circuit 12, an MRI pulse excitation system circuit 14 and an MRE signal measurement 
system circuit 16. The MRI control system circuit 12 controls operations of the fMRI 
system 10 as described herein. The MRI pulse excitation system circuit 14 and the ♦ 
MRI signal measurement system circuit 16 are controlled as described herein so as to 

20 provide a single excitation pulse and corresponding gradient and compensation pulses 
as well as provide for the partial £-space acquisition of the resulting MRI signals. 
Conventional MRI systems, such as those provided by General Electric Medical 
Systems, Siemens, Philips, Varian, Bruker, Marconi and Toshiba may be modified, by 
for example, programming such systems, to carry out the operations as described 

25 herein and, thereby, provide fMRI systems according to embodiments of the present 
invention. 

While an exemplary fMRI system is illustrated in Figure 1 and described 
herein with a particular division of functions and/or operations, as will be appreciated 
by those of skill in the art, other divisions of functions and/or operations may be 
30 utilized while still benefiting from the teachings of the present invention. For 

example, the MRI control system circuit 12 could be combined with either the MRI 
pulse excitation system circuit 14 or the MRI signal measurement system circuit 16. 
Thus, the present invention should not be construed as limited to a particular 
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architecture or division of MRI functions/operations but is intended to cover any 
architecture or division of functions/operations capable of carrying out the operations 
described herein. 

Operations according to embodiments of the present invention are illustrated 
5 in the flowchart of Figure 2 A, which will be described herein with reference to the 
exemplary fMRI system 10 of Figure 1. As seen in Figure 2A, the MRI control 
system circuit 12 causes the MRI pulse excitation system circuit 14 to excite a target 
region(s) of a subject with a single excitation pulse (block 30). The MRI control 
system circuit 12 also receives signal information from the MRI signal measurement 

10 system circuit 16 and obtains a first image using a partial A:-space frame from of the 
target based on the received signal information (block 32). After obtaining the first 
image, the MRI control system circuit 12 causes the MRI pulse excitation system 
circuit 14 to subject the target region(s) to a compensation pulse (block 34). The MRI 
control system circuit 12 receives signal information from the MRI signal 

15 measurement system circuit 16 subsequent to the target region(s) being subjected to 
the compensation pulse and obtains a second, compensated, image using a partial k- 
space frame from of the target based on the received signal information (block 36). 
The two images are then combined to provide the fMRI image (block 38). 

Figure 2B is a flowchart illustrating particular embodiments of the present 

20 invention. As seen in Figure 2B, axial slice selection can be provided by generating a 
matched quadratic excitation pulse (block 50) and exciting the target region with the 
matched quadratic excitation pulse (block 52). After axial slice selection using the 
matched quadratic excitation pulse, the first image is collected using a partial £-space 
acquisition in a first trajectory (e.g. a first half &-space plus additional seven more 

25 lines) (block 54). The target is then subjected to a brief instance of z-shimming 

gradient (block 56). A second image is collected using a second partial £-space frame 
(e.g. a second half &-space plus additional seven more lines) (block 58). The £-space 
trajectory of the second A:-space frame preferably has the same sampling direction to 
ensure the same characteristics of the phase propagation. Also, the center portions of 

30 the two partial A>space frames (kz frames) can, in certain embodiments, be kept 
substantially close together to reduce the T2* contrast mixing and also to provide 
sufficient signal level for the second acquisition. Finally, the first and second images 
are combined to provide the resultant image (block 60). 
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In certain embodiments of the present invention, an additional line can be 
positioned in the center (and/or at other desired locations) for both kz frames to serve 
as a reference to realign the even and odd lines in fc-space and correct for the phase 
errors. Such a realignment and phase correction may be carried out as described in 
5 E. C. Wong, "Shim-insensitive phase correction for EPI using a two-echo reference 
scan," Proc. SMRM, p. 4514, 1992. 

A schematic illustration of a pulse sequence according to embodiments of the i 
present invention is shown in Figure 3 A; the resultant £-space trajectory is shown in 
Figure 3B. As is seen in Figure 3A, the target area of a subject is subjected to a 

10 single matched quadratic excitation pulse 100. A slice is selected with the gradient 
Gz pulse 106 and a Gx gradient pulse 102 and a Gy gradient pulse 104 are also 
provided to as to establish the initial starting location of the initial partial &-space 
acquisition. The two kz frames are offset by Afe which is determined by the z- 
shimming gradient 110 and the initial starting location for the second partial k-space"^ 

15 acquisition is established by the Gy gradient pulse 108. The resulting trajectory is *' 
seen in Figure 3B. 

As is seen in Figure 3A, as a result of the sampling of the first and second 
partial Ar-space frames, the centers of the fc-space frames are corresponding close to the 
end of the sampling of the first partial £-space frame and close to the beginning of the 

20 sampling of the second partial £-space frame. As a result, Echo Time 1, the time 
difference between the excitation pulse and the center of the first partial k-space, 
(TEi) for the first partial k-space frame and Echo Time 2, the time difference between 
the excitation pulse and the center of the second partial k-space) (TE2) for the second 
partial k-space frame are relatively close together. Thus, according to embodiments 

25 of the present invention, it is preferred that TEi and TE2 be asymmetric and mirrored 
about the division between the first and second partial £-space frames. Furthermore, 
it is preferred that TEi and TE2 be close together to minimize signal decay from the 
first partial k-space acquisition to the second, preferably within about 10 ms. 

In Figure 3B, the first half of the first k frame 120 (shown as a solid line) and 

30 the second half of the second kz frame 124 (shown as a dashed line) are completely 
sampled, while the second half of the first kz frame 120 and first half of the second kz 
frame 124 are partially sampled to help arrange the images such that the centers of the 
£-space of both kz planes remain proximate and close together. The first kz frame 120 



11 



WO 02/35250 



r CT7US01/32277 



has an initial trajectory in the direction shown by arrow 122. The second kz frame 124 
preferably also has an initial trajectory in the same direction as the arrow 122 as is 
illustrated by the arrow 126. Thus, as described above, the first kz frame 120 and the 
second kz frame 124 may provide the same direction of sampling. Finally, the 
5 reference center lines are illustrated as lines 121 and 125 in Figure 3B. 

As will be appreciated by those of skill in the art in light of the present 
disclosure, partial k-space acquisition of a single excitation pulse may be provided 
through programming of an MRI system. Such programming may take many forms 
depending on the particular system utilized. However, set forth below is an example 
10 of a generalized program which provides a single excitation pulse and corresponding 
partial fc-space acquisition according to particular embodiments of the present 
invention. 

EXTERNAL_PULSE(rf1 , magnitude_rf1, rf Imposition, rf1_duration, 
slicethickness, flip_angle, pulse_cycle, number_of_points, 
1 5 quadratic_wave_magnitude); 

EXTERNAL_PULSE(phase, phase_rfl, rf1_position, rf1_duration, 
slicethickness, flip_angle, pulse__cycle, number_of_points, 
quadratic_wave_phase); 

20 

TRAPEZOID(GRADlENT_Z, gz1 jrephase, gz1_position, gz1_duration, 
gz1_amplitude); 

TRAPEZOID(GRADlENT_X, gxep_dephase, dephase_position, 
25 dephase_duration, gxep_dephase_amp1itude); 

for (i=0; i<epyres; { 

if (i%2 == 0) TRAPEZOID(GRADIENT_X, gxep_positive, 
start_position + i*cyclejength, gxep_duration, gxep_amplitude); 
30 else TRAPEZOID(GRADIENT_X, gxep_negative, start_position + 

i*cycle_Iength, gxep__duration, gxep_amplitude_negative); 
} 

TRAPEZOID(GRADIENT_Y, gyep_dephase, dephase_j>osition, 
35 dephase_duration, gyep_dephase_amplitude); 

for (i=0, j=0; i<epyres-3; j++) { 
if (i<=36) { 
if (j==32) j+=1; 

40 TRIANGLE(GRADIENT_Y, gyep_blip, blip_position + fcyclejength, 

blip_duration, blip_amplitude); 

if (i>37) { 
if(i==43)j+=1; 

12 
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TRIANGLE(GRADIENT_Y, gyep_blip, blip_position + j*cyclejength, 
blip_du ration, blip_amplitude); 
} 

5 TRAPEZOID(GRADIENT_Z, gzblip, position_gzblip, duration_gzblip, 

gzblip__amplitude); 

TRAPEZOID(GRADIENT_Y, gyblip, position_jgyblip, duration_gyb!ip, 
gyblip_amplitude); 

1 0 The above code has been generalized differing routines may be utilized 

depending on the particular MRI system utilized. However, in general, i the above 
code, EXTERN ALJPULSE is a routine that displays a sequence of binary numbers 
according to the arguments such as timing, magnitude and phase information of a 
given file, e.g., quadratic_wave_magnitude. TRAPEZOID is a routine that displays a 

15 trapezoid waveform according to the location, duration and amplitude given in the 
arguments. TRIANGLE is a routine that displays a triangle waveform according to 
the location, duration and amplitude given in the arguments. The above program is ~_ 
provided merely as an example. The present invention should not be construed as 
limited to such a program. 

20 Embodiments of the present invention may provide an effective single-shot' 

method that can recover susceptibility effect-induced signal losses. Compared with- 
previous methods that used multiple excitations, embodiments of the present 
invention can achieve comparable signal recovery with much improved temporal 
resolution. This advantage can allow fMRI investigations to be carried out in 

25 statically inhomogeneous areas such as the ventral frontal and inferior temporal areas 
without compromising acquisition speed. In certain embodiments and applications, it 
is anticipated that this technique can be particular useful in event-related studies of the 
brain, including human memory and attention processes that involve areas near 
air/tissue interfaces. Examples of such processes may include the progression of 

30 Alzheimer's disease, stroke, age, attention deficit disorders, and the like. In addition 
the impact of pharmacological products on this region of the brain may be verified or 
studied (before and after, or during a treatment) to ascertain the impact of a course of 
treatment or a targeted drug. 

Embodiments of the present invention will now be illustrated with reference to 

35 certain examples which are included herein for the purposes of illustration only, and 
which are not intended to be limiting of the invention. 
Examples 
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Imaging parameters 

A pulse sequence according to embodiments of the present invention was 
developed and implemented on a GE 1 .5 T whole-body scanner (General Electric 
Medical Systems, Milwaukee, Wisconsin). A complex excitation pulse with 
5 quadratic phase profile across the slice was used to better match the susceptibility- 
induced gradient. Raw data before Fourier Transform was collected using the 
aforementioned £-space trajectory. The slice thickness was 5 mm and FOV 24 cm. 
Both kz frames were acquired using a 64 x 39 matrix. The repetition time (TR) was 1 
second. The joint point of the two frames occurred at 40 ms from the excitation, 

1 0 resulting in echo times (TE) of 36 ms and 44 ms for the two images, respectively. A 
receiver bandwidth of 125 kHz was used, the transition time from zero to readout 
amplitude for the readout gradient was 96 \is. Thus, the total acquisition time was 55 
ms, comparable to a 45 ms acquisition window for a conventional single-shot 
gradient-recalled EPI of 64 x 64 matrix using the same slew rate. The z-shim gradient 

1 5 had a duration of 1 .4 ms at the joint of the two kz frames. Its amplitude could be 

altered by a control variable in real time to ensure optimal compensation. With this 
setup, a maximum of 13 slices can be acquired per second, highly comparable to a 
conventional single-shot EPI scan of 14 slices per second. 



raw data was acquired and transferred to a Silicon Graphics SGI 02 workstation 
(Mountain View, CA) for reconstruction off-line. It was first split into two parts from 
the joint of the two kz frames. The two partial kz frames were then squared by 
conjugate mirroring to reduce the discontinuity in the point-spread function. Through 

25 a realignment and phase correction procedure via the center reference lines, images 
were reconstructed by a conventional fast Fourier Transform. With z-shim gradient 
turned off, axial images of a water phantom from these two partial k-space 
acquisitions were obtained at a TR of 1 s and TE's of 36 ms and 44 ms. Figures 4A 
and 4B show all eight slices with slice thickness of 5 mm from these two acquisitions. 

30 The signal-to-noise ration (SNR) of the second acquisition was 12% less than that of 
the first acquisition, caused by the T2* decay. The second acquisition also showed 
slightly more distortion caused by the propagation of the accumulative phase errors. 
The resultant two images can then be summed to a final image by a simple magnitude 



Results 



20 



Experiments were carried out on a GE 1 .5 T whole-body MRI scanner. The 
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addition. To assess the feasibility of the in vivo application of this technique, images 
were also acquired in human subjects. Eight slices of axial images were acquired at 
the same TR of 1 s and TEs of 36 ms and 44 ms. The slice thickness was 5 mm. After 
a calibration run where a collection of ten z-shim gradients (thus, the shapes of the 
5 quadratic profile) were used, an optimal z-shim gradient was determined at 0.15 
G/cm. This procedure takes a period of ten TRs and only consumes short amount of 
time. At the optimal gradient amplitude, combined images of the two acquisitions are 
shown in Figure 5C. To serve as comparisons, images from the same locations using 
a conventional single-shot gradient-recalled EPI method at a TE of 40 ms are shown 

10 in Figure 5 A, and combined images using a two-shot gradient-recalled EPI method at 
a TE of 40 ms and with the same z-shim gradient are shown in Figure 5B. The 
average SNR of Figure 5B and Figure 5C are highly comparable at 103 and 101 
respectively; however, both are higher than that of 92 in Figure 5A. The reasons that 
the partial £-space acquisitions did not suffer much SNR may be mainly the result of ^ . 

15 completely sampled center portions of the two partial Ar-space, and also the slightly 
shorter TE for the first partial fc-space acquisition. Several areas at the ventral brain 
region are shown to have recovered signals (indicated by the arrows in Figures 5A 
through 5C). It can be seen that relatively uniform images are acquired without 
sacrificing temporal resolution. 

20 The present invention finds use for both pre-clinical animal studies, veterinary 

and medical applications. The present invention may be advantageously employed 
for diagnostic evaluation and/or treatment of subjects, in particular human subjects. 
"Subjects" according to the present invention can be any animal subject, and are 
preferably mammalian subjects (e.g., humans, canines, felines, bovines, caprines, 

25 ovines, equines, rodents, porcines, and/or lagomorphs), and more preferably are 
human subjects. 

The present invention may further be employed for: assessment of cerebral 
perfusion in following induced subarachnoid hemorrhage or in conditions marked by 
brain dysfunction, e.g., in connection with acute severe symptomatic hyponatremia; 
30 evaluation of new therapies, e.g. , in the treatment of cerebral vasospasm (including 
but not limited to, anti-thrombolytic therapies, calcium channel blockers, anti- 
inflammatory therapies, angioplasty, and the like); assessment of the presence or 
absence and/or severity of ischemia in large tissue masses; assessment of the 
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relationship between blood metabolites and cerebral perfusion in cerebral ischemia 
associated with acute liver failure, e.g., for the treatment of Alzheimer's disease; 
evaluation of new therapies for stroke, including but not limited to, t-PA, aspirin 
antiphospholipids, lupus anticoagulants, antiphospholipid antibodies, and the like; 
5 evaluation of risk factors for stroke, e.g. , serum lipid levels; evaluation of induced 
brain hypothermia on cerebral perfusion during neurosurgery, e.g., for stroke; 
evaluation of the effects of age on cerebral perfusion, e.g., to study lacunar infarcts; 
and assessment of narcotics, e.g., cocaine, amphetamines, ethanol, and the like, on the 
ischemic brain. 

10 The foregoing is illustrative of the present invention and is not to be construed 

as limiting thereof. Although a few exemplary embodiments of this invention have 
been described, those skilled in the art will readily appreciate that many modifications 
are possible in the exemplary embodiments without materially departing from the 
novel teachings and advantages of this invention. Accordingly, all such modifications 

15 are intended to be included within the scope of this invention as defined in the claims. 
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THAT WHICH IS CLAIMED IS: 

1 . A method of recovering a signal of a functional magnetic resonance 
imaging system, comprising: 

exciting a target region of a subject with a single excitation pulse; then 
obtaining a first image using a first partial k-space frame of the target region; 

then 

subjecting the target region to a compensation pulse; then 
obtaining a second, compensated, image using a second partial k-space frame 
of the target region; and 

combining the first and second images to form a combined image of the target 

region. 

2. The method of Claim 1, wherein the single excitation pulse comprises 
a matched quadratic excitation pulse. 

3. The method of Claim 2; wherein the compensation pulse comprises a 
z-shimming pulse. 

4. The method of Claim 1, wherein the first partial k-space frame has 
associated therewith a first sampling direction and the second partial k-space frame 
has associated therewith a second sampling direction and wherein the first sampling 
direction and the second sampling direction are substantially the same direction. 

5. The method of Claim 1 , wherein the first partial k-space frame and the 
25 second partial k-space frame are sampled so that a center of each partial k-space 

frame are proximate each other. 

6. The method of Claim 1 , wherein the step of obtaining a first image 
comprises the steps of: 

30 completely sampling a first half of the first partial k-space frame; and 

partially sampling a second half of the first partial k-space frame; and 
wherein the step of obtaining a second image comprises the steps of: 
partially sampling a first half of the second k-space frame; and 



5 



10 



15 



20 



17 



WO 02/35250 



• 




r CT/US01/32277 



completely sampling a second half of the second k-space frame. 



7. 



The method of Claim 6, wherein sampling of the first half and the 



second half of the first k-space frame and sampling of the first half and the second 
5 half of the second k-space frame are carried out following substantially the same 



8. The method of Claim 1 , wherein a TE associated with the first partial 
k-space frame and a TE associated with the second partial k-space frame are 
10 asymmetric with reference to the beginning of the respective k-space frames and are 
mirrored about a division between the first and second partial k-space frames. 



alignment indicia in the first and second images, and wherein the combining step 
1 5 comprises the step of aligning the indicia in each of the first and second images. 

10. The method of Claim 9, wherein the alignment indicia comprise a 
centerline of the first k-space frame and a centerline of the second k-space frame. 



sampling a centerline of the first partial k-space frame in two sampling 
directions, wherein the centerline comprises the center of a full k-space frame 
corresponding to the first partial k-space frame; and 

sampling a centerline of the second partial k-space frame in two sampling 
25 directions, wherein the centerline comprises the center of a full k-space frame 
corresponding to the second partial k-space frame. 

12. The method of Claim 1 , wherein the target region of a subject is an in 
vivo ventral frontal or inferior temporal area of the human brain. 



direction. 



9. 



The method of Claim 1, further comprising the step of inserting 



20 



11. 



The method of Claim 1, further comprising the steps of: 



30 



13 . The method of Claim 1 , wherein said combined image provides 
information about the function of the human brain including human memory and 
attention processes in areas near air/tissue interfaces. 
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14. A system for recovering a signal of a functional magnetic resonance 
imaging system, comprising: 

means for exciting a target region of a subject with a single excitation pulse; 
5 means for obtaining a first image using a first partial k-space frame of the 

target region; 

means for subjecting the target region to a compensation pulse responsive to 
the means for obtaining a first image; 

means for obtaining a second, compensated, image using a second partial k- 
10 space frame of the target region responsive to the means for exciting the target region 
with a compensation pulse; and 

means for combining the first and second images to form a combined image of 
the target region. 

15 15. The system of Claim 14, wherein the single excitation pulse comprises 

a matched quadratic excitation pulse. 

1 6. The system of Claim 15, wherein the compensation pulse comprises a 
z-shimming pulse. 

20 

17. The system of Claim 14, wherein the first partial k-space frame has 
associated therewith a first sampling direction and the second partial k-space frame 
has associated therewith a second sampling direction and wherein the first sampling 
direction and the second sampling direction are substantially the same direction. 

25 

1 8. The system of Claim 14, wherein the first partial k-space frame and the 
second partial k-space frame are sampled so that a center of each partial k-space 
frame is proximate to each other. 

30 19. The method of Claim 14, wherein the means for obtaining a first image 

comprises: 

means for completely sampling a first half of the first partial k-space frame; 

and 
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means for partially sampling a second half of the first partial k-space frame; 

and 

wherein the means for obtaining a second image comprises: 

means for partially sampling a first half of the second k-space frame; and 



20. The system of Claim 19, wherein sampling of the first half and the 
second half of the first k-space frame and sampling of the first half and the second 
half of the second k-space frame are carried following substantially the same 

10 direction. 

21 . The system of Claim 14, wherein a TE associated with the first partial 
k-space frame and a TE associated with the second partial k-space frame are 
asymmetric with reference to the beginning of the respective k-space frames and are 

1 5 mirrored about a division between the first and second partial k-space frames. 

22. The system of Claim 14, further comprising means for inserting 
alignment indicia in the first and second images, and wherein the means for 
combining comprises means for aligning the indicia in each of the first and second 

20 images. 

23. The system of Claim 22, wherein the alignment indicia comprise a 
centerline of the first k-space frame and a centerline of the second k-space frame. 

25 24. The system of Claim 14, further comprising: 

means for sampling a centerline of the first partial k-space frame in two 
sampling directions, wherein the centerline comprises the center of a full k-space 
frame corresponding to the first partial k-space frame; and 

means for sampling a centerline of the second partial k-space frame in two 
30 sampling directions, wherein the centerline comprises the center of a full k-space 
frame corresponding to the second partial k-space frame. 



5 



means for completely sampling a second half of the second k-space frame. 
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25. The system of Claim 14, wherein the target region of a subject is an in 
vivo ventral frontal or inferior temporal area of the human brain. 

26. The system of Claim 14, wherein the combined image provides 

5 information about the function of the human brain including human memory and 
attention processes in areas near air/tissue interfaces. 

27. A computer program product for recovering a signal of a functional 
magnetic resonance imaging system, comprising: 

10 a computer readable storage medium having computer readable program code 

embodied therein, the computer readable program code comprising: 

computer readable program code configured to excite a target region of a 
subject with a single excitation pulse; 

computer readable program code configured to obtain a first image using a - :r 
15 first partial k-space frame of the target region; 

computer readable program code configured to subject the target region to a 
compensation pulse responsive to the means for obtaining a first image; 

computer readable program code configured to obtain a second, compensated, 
image using a second partial k-space frame of the target region responsive to the 
20 means for exciting the target region with a compensation pulse; and 

computer readable program code configured to combine the first and second 
images to form a combined image of the target region. 

28. The computer program product of Claim 27, wherein the single 
25 excitation pulse comprises a matched quadratic excitation pulse. 

29. The computer program product of Claim 28, wherein the compensation 
pulse comprises a z-shimming pulse. 

30 30. The computer program product of Claim 27, wherein the first partial k- 

space frame has associated therewith a first sampling direction and the second partial 
k-space frame has associated therewith a second sampling direction and wherein the 
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first sampling direction and the second sampling direction are substantially the same 
direction. 



3 1 . The computer program product of Claim 27, wherein the first partial k- 
5 space frame and the second partial k-space frame are sampled so that a center of each 

partial k-space frame is proximate each other. 

32. The computer program product of Claim 27, wherein the means for 
obtaining a first image comprises: 

computer readable program code configured to completely sample a first half 
of the first partial k-space frame; and 

computer readable program code configured to partially sample a second half 
of the first partial k-space frame; and 

wherein the means for obtaining a second image comprises: 
computer readable program code configured to partially sample a first half of 
the second k-space frame; and 

computer readable program code configured to completely sample a second 
half of the second k-space frame. ' 

20 33. The computer program product of Claim 32, wherein sampling of the 

first half and the second half of the first k-space frame and sampling of the first half 
and the second half of the second k-space frame are carried out in substantially a same 
direction. 



10 



15 



25 34. The computer program product of Claim 27, wherein a TE associated 

with the first partial k-space frame and a TE associated with the second partial k- 
space frame are asymmetric with reference to the beginning of the respective k-space 
frames and are mirrored about a division between the first and second partial k-space 
frames. 

30 

3 5 . The computer program product of Claim 27, further comprising 
computer readable program code which inserts alignment indicia in the first and 
second images, and wherein the computer readable program code configured to 
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combine comprises computer readable program code configured to align the indicia in 
each of the first and second images. 

36. The computer program product of Claim 35, wherein the alignment 

, 5 indicia comprise a centerline of the first k-space frame and a centerline of the second 
k-space frame. 

37. The computer program product of Claim 27, further comprising: 
computer readable program code configured to sample a centerline of the first 

10 partial k-space frame in two sampling directions, wherein the centerline comprises the 
center of a full k-space frame corresponding to the first partial k-space frame; and 
computer readable program code configured to sample a centerline of the 
second partial k-space frame in two sampling directions, wherein the centerline 
comprises the center of a full k-space frame corresponding to the second partial k- «V 

15 space frame. 

38. The computer program product of Claim 27, wherein the target region 
of a subject is an in vivo ventral frontal or inferior temporal area of the human brain. 

20 39. The computer program product of Claim 27, wherein the combined 

image provides information about the function of the human brain including human 
memory and attention processes in areas near air/tissue interfaces. 
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